Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria).

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

COI metabarcoding of large benthic Foraminifera: Method validation for application in ecological studies.

Monitoring community composition of Foraminifera (single-celled marine protists) provides valuable insights into environmental conditions in marine ecosystems. Despite the efficiency of environmental DNA (eDNA) and bulk-sample DNA (bulk-DNA) metabarcoding to assess the presence of multiple taxa, this has not been straightforward for Foraminifera partially due to the high genetic variability in widely used ribosomal markers. Here, we test the correctness in retrieving foraminiferal communities by metabarcoding of mock communities, bulk-DNA from coral reef sediment samples, and eDNA from their associated ethanol preservative using the recently sequenced cytochrome c oxidase subunit 1 (COI) marker. To assess the detection success, we compared our results with large benthic foraminiferal communities previously reported from the same sampling sites. Results from our mock communities demonstrate that all species were detected in two mock communities and all but one in the remaining four. Technical replicates were highly similar in number of reads for each assigned ASV in both the mock communities and bulk-DNA samples. Bulk-DNA showed a significantly higher species richness than their associated eDNA samples, and also detected additional species to what was already reported at the specific sites. Our study confirms that metabarcoding using the foraminiferal COI marker adequately retrieves the diversity and community composition of both the mock communities and the bulk-DNA samples. With its decreased variability compared with the commonly used nuclear 18 S rRNA, the COI marker renders bulk-DNA metabarcoding a powerful tool to assess foraminiferal community composition under the condition that the reference database is adequate to the target taxa.

Mitochondrial cytochrome c oxidase subunit I (COI) metabarcoding of Foraminifera communities using taxon-specific primers.

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in most marine environments. Molecular methods such as metabarcoding have revealed a high, yet undescribed diversity of Foraminifera. However, so far only one molecular marker, the 18S ribosomal RNA, was available for metabarcoding studies on Foraminifera. Primers that allow amplification of foraminiferal mitochondrial cytochrome oxidase I (COI) and identification of Foraminifera species were recently published. Here we test the performance of these primers for the amplification of whole foraminiferal communities, and compare their performance to that of the highly degenerate LerayXT primers, which amplify the same COI region in a wide range of eukaryotes. We applied metabarcoding to 48 samples taken along three transects spanning a North Sea beach in the Netherlands from dunes to the low tide level, and analysed both sediment samples and meiofauna samples, which contained taxa between 42 µm and 1 mm in body size obtained by decantation from sand samples. We used single-cell metabarcoding (Girard et al., 2022) to generate a COI reference library containing 32 species of Foraminifera, and used this to taxonomically annotate our community metabarcoding data. Our analyses show that the highly degenerate LerayXT primers do not amplify Foraminifera, while the Foraminifera primers are highly Foraminifera- specific, with about 90% of reads assigned to Foraminifera and amplifying taxa from all major groups, i.e., monothalamids, Globothalamea, and Tubothalamea. We identified 176 Foraminifera ASVs and found a change in Foraminifera community composition along the beach transects from high tide to low tide level, and a dominance of single-chambered monothalamid Foraminifera. Our results highlight that COI metabarcoding can be a powerful tool for assessing Foraminiferal communities.

First report of mitochondrial COI in foraminifera and implications for DNA barcoding.

Foraminifera are a species-rich phylum of rhizarian protists that are highly abundant in many marine environments and play a major role in global carbon cycling. Species recognition in Foraminifera is mainly based on morphological characters and nuclear 18S ribosomal RNA barcoding. The 18S rRNA contains variable sequence regions that allow for the identification of most foraminiferal species. Still, some species show limited variability, while others contain high levels of intragenomic polymorphisms, thereby complicating species identification. The use of additional, easily obtainable molecular markers other than 18S rRNA will enable more detailed investigation of evolutionary history, population genetics and speciation in Foraminifera. Here we present the first mitochondrial cytochrome c oxidase subunit 1 (COI) gene sequences ("barcodes") of Foraminifera. We applied shotgun sequencing to single foraminiferal specimens, assembled COI, and developed primers that allow amplification of COI in a wide range of foraminiferal species. We obtained COI sequences of 49 specimens from 17 species from the orders Rotaliida and Miliolida. Phylogenetic analysis showed that the COI tree is largely congruent with previously published 18S rRNA phylogenies. Furthermore, species delimitation with ASAP and ABGD algorithms showed that foraminiferal species can be identified based on COI barcodes.

The mitochondrial genome of Nemalecium lighti (Hydrozoa, Leptothecata).

The hydrozoan species Nemalecium lighti (Hargitt, 1924) is widely distributed in tropical marine waters around the world. Here we report the complete linear mitochondrial genome of N. lighti from Sint Eustatius (Lesser Antilles). The mitochondrial genome with a length of 14,320 bp encodes for 13 protein-coding genes, two tRNA genes, and two rRNA genes. Gene arrangement differs from that found in other species of the same taxonomic order and a phylogenetic analysis shows that based on mitochondrial genes, N. lighti clusters outside of the Leptothecata, rendering the order paraphyletic.

Integrating morphology and metagenomics to understand taxonomic variability of Amphisorus (Foraminifera, Miliolida) from Western Australia and Indonesia.

Foraminifera are a group of mostly marine protists with high taxonomic diversity. Species identification is often complex, as both morphological and molecular approaches can be challenging due to a lack of unique characters and reference sequences. An integrative approach combining state of the art morphological and molecular tools is therefore promising. In this study, we analysed large benthic Foraminifera of the genus Amphisorus from Western Australia and Indonesia. Based on previous findings on high morphological variability observed in the Soritidae and the discontinuous distribution of Amphisorus along the coast of western Australia, we expected to find multiple morphologically and genetically unique Amphisorus types. In order to gain detailed insights into the diversity of Amphisorus, we applied micro CT scanning and shotgun metagenomic sequencing. We identified four distinct morphotypes of Amphisorus, two each in Australia and Indonesia, and showed that each morphotype is a distinct genotype. Furthermore, metagenomics revealed the presence of three dinoflagellate symbiont clades. The most common symbiont was Fugacium Fr5, and we could show that its genotypes were mostly specific to Amphisorus morphotypes. Finally, we assembled the microbial taxa associated with the two Western Australian morphotypes, and analysed their microbial community composition. Even though each Amphisorus morphotype harboured distinct bacterial communities, sampling location had a stronger influence on bacterial community composition, and we infer that the prokaryotic community is primarily shaped by the microhabitat rather than host identity. The integrated approach combining analyses of host morphology and genetics, dinoflagellate symbionts, and associated microbes leads to the conclusion that we identified distinct, yet undescribed taxa of Amphisorus. We argue that the combination of morphological and molecular methods provides unprecedented insights into the diversity of foraminifera, which paves the way for a deeper understanding of their biodiversity, and facilitates future taxonomic and ecological work.

Uncovering bacterial and functional diversity in macroinvertebrate mitochondrial-metagenomic datasets by differential centrifugation.

PCR-free techniques such as meta-mitogenomics (MMG) can recover taxonomic composition of macroinvertebrate communities, but suffer from low efficiency, as >90% of sequencing data is mostly uninformative due to the great abundance of nuclear DNA that cannot be identified with current reference databases. Current MMG studies do not routinely check data for information on macroinvertebrate-associated bacteria and gene functions. However, this could greatly increase the efficiency of MMG studies by revealing yet overlooked diversity within ecosystems and making currently unused data available for ecological studies. By analysing six 'mock' communities, each containing three macroinvertebrate taxa, we tested whether this additional data on bacterial taxa and functional potential of communities can be extracted from MMG datasets. Further, we tested whether differential centrifugation, which is known to greatly increase efficiency of macroinvertebrate MMG studies by enriching for mitochondria, impacts on the inferred bacterial community composition. Our results show that macroinvertebrate MMG datasets contain a high number of mostly endosymbiont bacterial taxa and associated gene functions. Centrifugation reduced both the absolute and relative abundance of highly abundant Gammaproteobacteria, thereby facilitating detection of rare taxa and functions. When analysing both taxa and gene functions, the number of features obtained from the MMG dataset increased 31-fold ('enriched') respectively 234-fold ('not enriched'). We conclude that analysing MMG datasets for bacteria and gene functions greatly increases the amount of information available and facilitates the use of shotgun metagenomic techniques for future studies on biodiversity.

Assessing the influence of sample tagging and library preparation on DNA metabarcoding.

Metabarcoding is increasingly used to assess species diversity by high-throughput sequencing where millions of sequences can be generated in parallel and multiple samples can be analysed in one sequencing run. Generating amplified fragments with a unique sequence identifier ('tag') for each sample is crucial, as it allows assigning sequences to the original samples. The tagging through so-called fusion primers is a fast and cheap alternative to commercially produced ligation-based kits. However, little is known about potential bias and inconsistencies introduced by the long nucleotide tail attached to those primers, which could lead to deficient reports of community composition in metabarcoding studies. We therefore tested the consistency and taxa detection efficiency of fusion primers in (1) a one-step and (2) two-step PCR protocol as well as (3) a commercially manufactured Illumina kit using mock communities of known composition. The Illumina kit delivered the most consistent results and detected the highest number of taxa. However, success of the two-step PCR approach was only marginally lower compared to the kit with the additional advantage of a much more competitive price per library. While most taxa were also detected with the one-step PCR approach, the consistency between replicates including read abundance was substantially lower. Our results highlight that method choice depends on the precision needed for analysis as well as on economic considerations and recommend the Illumina kit to obtain most accurate results and the two-step PCR approach as a much cheaper yet very robust alternative.

Comparison of environmental DNA and bulk-sample metabarcoding using highly degenerate cytochrome c oxidase I primers.

Freshwater biodiversity provides important ecosystem services and is at the core of water quality monitoring worldwide. To assess freshwater biodiversity, genetic methods such as metabarcoding are increasingly used as they are faster and allow better taxonomic resolution than manual identification methods. Either sampled organisms are used directly for "bulk metabarcoding," or water is filtered and the extracted environmental DNA serves as a proxy for biodiversity via "eDNA metabarcoding." Despite the advantages of both methods, questions remain regarding their comparability and applicability for routine biomonitoring and stressor impact assessment. Therefore, we compared metabarcoding results from bulk and eDNA samples taken from 19 streams spanning a wide gradient of farming intensities in New Zealand. We performed PCR with highly degenerate cytochrome c oxidase I primers and sequenced libraries on an Illumina MiSeq. The inferred community composition differed strongly between the two methods. More taxa were captured by eDNA than bulk-sample metabarcoding (5,819 vs. 1,483), but more of the commonly used invertebrate bioindicator taxa (mayflies, stoneflies and caddisflies) were found in bulk (47) than eDNA samples (37). Catchment-wide and local land use impacts on communities were detected better by eDNA metabarcoding, especially for non-metazoan taxa. Our findings imply that bulk-sample metabarcoding resembles classical freshwater biomonitoring approaches better, as more indicator macroinvertebrate taxa are captured. However, eDNA metabarcoding might be better suited to infer the impact of stressors on stream ecosystems at larger scales, as many new and potentially more informative taxa are registered. We therefore suggest exploring both methods in future assessments of stream biodiversity.

The complete mitochondrial genome of a cryptic amphipod species from the Gammarus fossarum complex.

The freshwater amphipod Gammarus fossarum is widely distributed throughout Europe and an important species for stream biomonitoring. It is known to consist of several cryptic species. We here report the complete mitochondrial genome of G. fossarum clade 11/type B with a length of 15,989 bp, encoding for 13 protein-coding genes, 22 tRNA genes, and 2 rRNA genes. Protein-coding and ribosomal genes have a similar arrangement as in other gammarid amphipods. A phylogenetic analysis clarifies the placement of G. fossarum within the Gammaridae.

Assessing the phylogeographic history of the montane caddisfly Thremma gallicum using mitochondrial and restriction-site-associated DNA (RAD) markers.

Repeated Quaternary glaciations have significantly shaped the present distribution and diversity of several European species in aquatic and terrestrial habitats. To study the phylogeography of freshwater invertebrates, patterns of intraspecific variation have been examined primarily using mitochondrial DNA markers that may yield results unrepresentative of the true species history. Here, population genetic parameters were inferred for a montane aquatic caddisfly, Thremma gallicum, by sequencing a 658-bp fragment of the mitochondrial CO1 gene, and 12,514 nuclear RAD loci. T. gallicum has a highly disjunct distribution in southern and central Europe, with known populations in the Cantabrian Mountains, Pyrenees, Massif Central, and Black Forest. Both datasets represented rangewide sampling of T. gallicum. For the CO1 dataset, this included 352 specimens from 26 populations, and for the RAD dataset, 17 specimens from eight populations. We tested 20 competing phylogeographic scenarios using approximate Bayesian computation (ABC) and estimated genetic diversity patterns. Support for phylogeographic scenarios and diversity estimates differed between datasets with the RAD data favouring a southern origin of extant populations and indicating the Cantabrian Mountains and Massif Central populations to represent highly diverse populations as compared with the Pyrenees and Black Forest populations. The CO1 data supported a vicariance scenario (north-south) and yielded inconsistent diversity estimates. Permutation tests suggest that a few hundred polymorphic RAD SNPs are necessary for reliable parameter estimates. Our results highlight the potential of RAD and ABC-based hypothesis testing to complement phylogeographic studies on non-model species.